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Country’s economic, social and cultural development is 
mainly dependent on the performance of its highway 
structure. Selection of appropriate pavement type and related 
design method are vital for the improvement of pavement 
performance and its service life, and reduction in the initial 
and maintenance cost. The rigid pavement exposed to many 
distresses during its service life resulted due to the variation 
of traffic loading, material properties, and _ climatic 
conditions. The main objective of this project is to make a 
comparison between manual and computer design for rigid 
pavement structure under different loading, material 
properties, and temperature regimes. For manual design and 
computer design, “Westergaard Method” and “KENPAVE 
software” were used respectively. The stress analysis results 
revealed that edge stresses are higher as compared with 
interior and corner location, and stresses estimated at all 
locations with Westergaard method are significantly lower 
than stresses estimated with KENPAVE software. Results of 
sensitivity analysis showed that change in pavement 
thickness, material properties, and wheel load has a 
significant impact on developed stresses at different slab 
locations. 
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1. Introduction 


Currently, many highway Projects are under construction in Oman, for Example, BATINAH 
EXPRESSWAY CONSTRUCTION PROJECT and BIDBID - SUR DUAL CARRIAGEWAY 
PROJECT. Generally, highway projects are high-cost projects and consume a large amount of 
country’s development budget. In addition, country’s economic, social and cultural development 
is mainly dependent on the performance of its highway structure. Selection of appropriate 
pavement type and related design method are vital for the improvement of pavement 
performance and its service life, and reduction in the initial and maintenance cost. Pavement is 
an artificial surface laid over the ground to simplify the travel. The pavement structure must be 
durable and serviceable in order to resist the imposed traffic loading. The top surface of rigid 
pavement structure usually consists of Portland cement concrete (PCC). PCC has a high modulus 
of elasticity and rigidity, and do not flex noticeably under the application of traffic loading [1]. 


The rigid pavement exposed to many types of distresses during its service life that includes 
warping cracks, joint faulting, and punch-out’s failure. The continuous variation in traffic 
loading, material properties, and climatic conditions mainly contributes towards such distresses 
and results in reduced pavement performance and its service life. To improve the performance of 
pavement, it is required to evaluate the behavior of pavement structure under varying loading, 
material properties, and climatic conditions. The main objective of this study is to make a 
comparison for stresses and deflections at various locations between Westergaard method and 
KENPAVE software. 


2. Literature review 


Concrete pavement is susceptible to the combination of thermal and traffic stress. When the load 
applied on the surface of the rigid pavement, it causes flexure of the slab that produces tensile 
and compressive stresses. There are three types of traffic loading in a concrete slab: loading at 
corner, edge loading, and interior loading. According to Westergaard equation, the highest stress 
level was caused by traffic induce stress is edge loading stress due to lower thickness of the 
grouted macadam layer, which is less than half of the concrete thickness [2]. 


Traffic load is one of the major factors that affect pavement performance. Equivalent single axle 
load (ESAL) provides the typical value of the traffic loading experienced by a pavement 
structure in its whole life [3]. Temperature stresses developed in cement concrete due to daily 
variation in slab temperature gradient along the slab thickness, and seasonal temperature 
variation due to the overall change in slab temperature. Khana et al. (2014) state that cement 
concrete assumed to be identical and have uniform elastic properties with subgrade. The main 
factor considered in the design of rigid pavements is the structural strength of the concrete; 
therefore, slight differences in subgrade strength have little influence upon the structural capacity 
of the pavement [4]. Pavements are designed to distribute traffic induced stresses to the 
subgrade. For this reason, subgrade condition is consider as a principal factor in selecting the 
pavement structure, and before the pavement is plotted, the quality of the subgrade soils must be 
evaluated to ensure that it has sufficient strength to carry the expected traffic loads during the 
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design life of the pavement. The pavement must also be plotted to limit the expansion and loss of 
density of the subgrade soil [5]. 


Arora (2003) has reported that the Westergaard’s analysis is used for the design of rigid 
pavements. The stresses in the concrete slab are determined using Westergaard’s theory. 
Westergaard considered the rigid pavement as a thin elastic plate resting on soil subgrade. The 
slab deflection depends on the stiffness of the subgrade and the flexural strength of the slab. 
Therefore, the stress-deformation characteristics of a rigid pavement depend on the relative 
stiffness of the slab and the subgrade. There are other methods used to calculate rigid pavement 
stresses other than using Westergaard method. Finite element method is suggested as an 
alternative method, which can provide optimum and economical design because of the procedure 
of FEM and calculates stresses at each node. This method can estimate accurately the 
performance of rigid pavement under the design load [6]. 


3. Research methodology 


Figure 1 shows the design framework of this research work. It starts with the selection of design 
input variables. These variables include traffic, subgrade soil properties, PCC properties and 
climatic data of the selected site. Using these input pavement structure is designed using 
Westergaard method and computer program. 


Selection of Input Parameters 
properties 


Poisson's 


Ratio for subgrade (K) 
Concrete 


Traffic loading 


Design by AASHTO and 
Westergaard Method 


Design by Computer Program 


Comparison of results 
Conclusions and suggestions 


Fig.1. Design framework. 


Initially, the pavement was design using AASHTO method and obtained results were used as a 
reference for other methods. 
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AASHTO Method: In this method, the empirical equation is used to related observed or 
measurable phenomena with the required outcomes. The equation-1 is widely used for the design 
of rigid pavement and has the following form. 


APSI 
log 9 — ‘ 075 _ 
45-15) | 99 -0132p,)x logy (st ¥c, X05 - 1.132) 


log oWig)= Zp XS, + 7.35 xlog yg(D+ 1)- 0.06 + 


215.63(7)} D°? - 


In above equation: 


Wis =predicted number of 80 KN (18,000 Ib) ESALs 


ZR = standard normal deviate 

So = combined standard error of the traffic prediction and performance prediction 
D = slab depth 

Pt = terminal serviceability index 

APSI = difference between the initial design serviceability index, po, and 


the design terminal serviceability index, p; 


Cd  =drainage coefficient 

J = load transfer coefficient 

E, = elastic modulus of PCC 

k = modulus of subgrade reaction 


ESAL: in pavement design, it is a common practice to convert the damage from wheel loads of 
various magnitudes and repetitions of mixed traffic to damage for an equivalent number of 
“standard single axle” loads. This conversion makes the design calculations simple. The most 
commonly used equivalent load in the U.S. is the 18,000 Ib (80 KN) [7]. The equation-2 is used 
to calculate ESAL for the design period at a growth rate of 4% [8]. 


ESALp = (AADT); x (Lp) X T24. X Dp x Ep x 365 (2) 
Where; 

AADT = average annual daily traffic 

Lr = Lane factor converts directional trucks to the design lane trucks 


Toa = Percent heavy trucks during a 24-hour period with six tires or more 
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Dr = directional distribution factor 
Er = Equivalency Factor, which is the damage caused by one average heavy 
truck measured in 18-KIP (80-kN) ESALs 


Westergaard Method: Westergaard assumed a pavement slab to be a thin plate resting on a 
special subgrade, which is considered elastic in a vertical direction only. The pressure- 
deformation characteristics of a rigid pavement depend upon the relative stiffness of the slab and 
the subgrade [9]. 


In this method, equation are available for three loading conditions, 1.e. interior, edge, and corner. 
In addition, curling stresses are calculated at edge and interior location only. These equations are 
given below for each location. In this method, the wheel load is assumed circular with radius ‘a’. 


Equation for interior stress (o;) due to circular load (radius = a): 


2 
ga na +05 pe a 
2h a SOR. 


(3) 
Equation for interior deflection (6;) due to circular load: 
P 1 a) 
6,=— 1+/ In| — |+ 7 -1.25 G 
8kl 20 2b l 
(4) 
Edge Stress (Semi Circular) Westergaard Equations: 
Equation for edge stress (o-) due to semi-circular load: 
3 
ph I Ni ay sds 22 
m(3+ “yh 100ka, 3 l 
(5) 
Equation for edge deflection (6,) due to semi-circular load: 
0.5 7] 
pe Ble EE) 903334017 0)| & 
(E hk) Ao (6) 


a, = av2 
Edge Stress (Circular) Westergaard Equations: 


Equation for edge stress (,) due to circular load: 


P : 4p I- 
ae 3(1+ Ll) ia Eh +184 2h EBs risis2y {S| 
m3+U)h’| 100ka’ 3 2 : 


(7) 
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Equation for edge deflection (6-) due to circular load: 


Fae a an ae 1-(o76+044(2) 


© (ERR) (8) 
Corner Stress Westergaard Equations: 
Equation for Corner stress (¢-) due to circular load (radius = a): 
0.6 
3P V2a 
Oc > 1 5 
he U 
(9) 


Equation for Corner deflection (6,) due to circular load: 


a= 2) 11-089{ 2 ‘)| 
(10) 


Where in equation 3 to 10: 


y = Euler's constant (0.577215) 
P=wheel load 


k = modulus of subgrade reaction 
/= Poisson ratio 

a = radius of circular load 

/ = length of slab 

h= slab thickness 

E = elastics modulus of pec 


Curling Stress Westergaard Equations: 


Exterior o,= eee (11) 
Interior o;= font Gi nee (12) 
2 I-p 
Where; 
C=1- POS ceo (tanA+tanh/ ) (13) 
sin 2A + sinh 2 2 
t= a (Edge) = eae (Interior) 


eJ8 eJ8 
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L,W =Length, and Width of PCC Slab 

C, C1, C2 =coefficients for interior stress for the desired direction 
AT = Temperature differential in degree F 

ar = Coefficient of thermal expansion 


KENPAVE SOFTWARE: Kenpave program is used to determine stresses, strains, and 
deflections at different points in the pavement structure subjected to traffic loading and taking 
into account the material properties. General information is assigned to the program like; the 
number of slab layers to be designed and the data of materials in each layer, e.g. thickness, 
poisons ratio, unit weight, and elastic modulus. In addition, type of foundation used whether it is 
liquid, solid, or Winkler foundation, number of periods per year and number of load groups. In 
addition, loading data is entered. Finally, the program is run to get the output (stress and 
deflection). 


4. Results and discussion 


Accumulative ESAL is calculated for arterial urban rigid pavement based on 8662 AADT, which 
projected for 30 years from 2016 until 2047. The opening year of the project is 2017. Since it is 
two-way traffic, the directional distribution factor (DF) taken as 0.5 and the percent of heavy 
trucks during a 24-hour period (T24) = 25. The lane factor=75 and the EF for arterial urban rigid 
pavement = 2.02. 


4.1. Design by AASHTO method 


The values of input variables were taken from AASHTO design guide, and other references and 
pavement thickness was calculated. The input data includes: reliability = 90%, combined 
standard Error (So) = 0.4, initial Serviceability Index (po) = 4.5, Terminal serviceability index (py 
= 3, Elastic modulus of PCC (Ec) = 4,000,000 psi (27586 MPa), Modulus of Rupture (S) = 700 
psi (4.826 MPa), Drainage coefficient (Cd) = 1, Load transfer coefficient (J) = 3.2, Modulus of 
subgrade reaction (k) from CBR% = 304.5 pci (82.76 KPa/mm), Standard normal deviate (Zr) = 
- 1.282, and APSI = 1.5. The calculated thickness using equation-lwas 300 mm and this 
thickness was used as a reference for analysis using Westergaard method and Computer Program. 


4.2. Stresses and deflections calculation by Westergaard method and KENPAVE 


Table 2 presents the summary of stresses and deflection of Westergaard method and KENPAVE 
analysis. For the base case, using Westergaard’s method Curling stresses are higher at the interior 
of the slab as compared with slab edge. Traffic-induced stresses are higher at the edge are higher 
in comparison to other locations. In addition, for KENPAVE analysis results show that edge 
stresses are higher than the interior and corner stresses that are similar to a manual method. 
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Table 1 
Calculation of ESAL. 

Year AADT ESAL Accumulative ESAL Dr Tx4 LF EF 
--- Veh/day --- --- --- % --- --- 
2016 8662 598762 --- 0.5 25 0.75 2.02 
2017 9009 622712 622712 0.5 25 0.75 2.02 
2018 9369 647621 1270333 0.5 25 0.75 2.02 
2019 9744 673526 1943859 0.5 25 0.75 2.02 
2020 10134 700467 2644325 0.5 25 0.75 2.02 
2021 10539 728485 3372811 0.5 25 0.75 2.02 
2022 10961 757625 4130436 0.5 25 0.75 2.02 
2023 11399 787930 4918365 0.5 25 0.75 2.02 
2024 11855 819447 5737812 0.5 25 0.75 2.02 
2025 12329 852225 6590037 0.5 25 0.75 2.02 
2026 12822 886314 7476351 0.5 25 0.75 2.02 
2027 13335 921766 8398117 0.5 25 0.75 2.02 
2028 13869 958637 9356754 0.5 25 0.75 2.02 
2029 14424 996982 10353737 0.5 25 0.75 2.02 
2030 15000 1036862 11390598 0.5 25 0.75 2.02 
2031 15601 1078336 12468935 0.5 25 0.75 2.02 
2032 16225 1121470 13590404 0.5 25 0.75 2.02 
2033 16874 1166328 14756733 0.5 25 0.75 2.02 
2034 17548 1212982 15969715 0.5 25 0.75 2.02 
2035 18250 1261501 17231215 0.5 25 0.75 2.02 
2036 18980 1311961 18543176 0.5 25 0.75 2.02 
2037 19740 1364439 19907616 0.5 25 0.75 2.02 
2038 20529 1419017 21326633 0.5 25 0.75 2.02 
2039 21350 1475778 22802410 0.5 25 0.75 2.02 
2040 22204 1534809 24337219 0.5 25 0.75 2.02 
2041 23093 1596201 25933420 0.5 25 0.75 2.02 
2042 24016 1660049 27593469 0.5 25 0.75 2.02 
2043 24977 1726451 29319920 0.5 25 0.75 2.02 
2044 25976 1795509 31115429 0.5 25 0.75 2.02 
2045 27015 1867329 32982759 0.5 25 0.75 2.02 
2046 28096 1942023 34924782 0.5 25 0.75 2.02 
2047 29219 2019704 36944485 0.5 25 0.75 2.02 
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Table 2 
Westergaard method and KENPAVE software results for Stresses and Deflections. 
Westergaard method KENPAVE software 
= Traffic . Traffic . 
Posit fL 
i Induced Cuclite Deflection Induced Cusine Deflection 
Stresses Stresses 
Stresses Gam Stresses Ga 
mm mm 
(MPa) (MPa) (MPa) (MPa) 
Interior 0.598 1.198 0.057 0.815 1.256 0.232 
Edge 1.139 1.081 0.214 1.462 1.170 0.386 
Corner 0.787 --- 0.361 0.873 --- 0.285 


Figure 2 shows the comparison between Westergaard method and KENPAVE software results for 
stresses at different locations. The estimated stresses with Westergaard method are lower than the 
computer method at all locations. Traffic-induced interior, edge, and corner stresses are lower by 
36.3%, 28.35, and 11%, respectively. The difference for curling stresses is less. These results 
implied that Westergaard method underestimates the stresses. This is due to its limitation in 
considering the subgrade supporting conditions for different locations. As Bartosova mentioned 
in his study that the highest load stress done by Westergaard method and finite element method 
occurs at longitudinal and transverse edges of the slab, and Westergaard method results in lower 
stresses in comparison to Finite element method [10]. 


1.6 
al 1.256 
1.198 

12 1.39 1.081 
= 1 
a. 
= 0.815 0.787 
= Os 
£ 0.598 
s 046 

o4 

0.2 

oO 
Interior Curling Stress Curling Stress 
(interior) (edge) 


MKENPAVE OO Westergaard 


Fig. 2. Comparison of stresses between Westergaard and KENPAVE method. 
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4.3. Sensitivity analysis for stresses 


The sensitivity analysis was conducted in order to assess the variation in stresses at a different 
level of material properties and traffic loading. This analysis was conducted using Westergaard's 
method. 


Stresses at Different Pavement Thickness: 


Stresses were calculated at different locations (interior, edge, and corner) under various slab 
thicknesses. Figure 3(a) shows that increase in slab thickness results a significant reduction in 
stresses at all locations. Figure 3(b) shows the relation between curling stresses and pavement 
slab thickness. Curling stresses are maximum at slab depth of 200 mm and stress decreases with 
the increase of slab depth. Curling stresses are resulted due to temperature difference at top and 
bottom of the slab. This temperature variation across slab thickness is usually non-linear, and this 
variation is more at top and little at the bottom of slab [11]. Thus the differential temperature 
change results reduction in curling stresses with the increase of slab depth. 


Stresses at Different PCC Modulus of Elasticity: 


The result of loading and curling stresses at different PCC modulus of elasticity with slab 
thickness = 12in (300mm) are presented in Figure 3(c). Results show that change in different 
stresses with an increase of PCC modulus of elasticity is linear. As modulus of elasticity 
increases stresses also increases and it has the impact that is more significant on traffic-induced 
stresses at the interior location and curling stresses at edge location. 


Stresses at Different PCC Thermal Coefficient of Expansion: 


Stresses were calculated under different PCC thermal coefficient of expansion. There was no 
significant change in loading stresses. Increase in thermal expansion coefficient increases the 
curling stresses within the pavement structure. The results are shown below in figure 3(d). 


Stresses at Different Slab Length: 


Figure 3(e) show that curling stresses increase with the increase of slab length. This is because 
slab geometric dimensions are one of governing parameters in resulting curling stresses along 
with material thermal properties. The temperature differential is also affected by the increase of 
slab length, which results in higher stresses. These results imply that in order to minimize the 
curling stresses appropriate length needed to be selected. 


Stresses at Different Wheel Loads: 


The stresses were calculated with respect to 7000, 9000, 12000 Ib wheel loads. Figure 3(f) shows 
that by increasing the wheel load the corner and edge loading stresses increase but the interior 
stresses remain constant. This may be due to different supporting conditions at the bottom of the 
slab at interior, edge, and corner. There is no significant effect on curling stresses, as they remain 
constant under different wheel loads. 
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Fig. 3. Results of sensitivity analysis for stresses at different locations. 


5. Conclusions 


e The maximum loading stress calculated is 1.139MPa, which occurs at the edge. 


e The maximum curling stress occurs at the interior, which is 1.198MPa. 


e For both Westergaard and sensitivity analysis, the critical loading stress appears at the 
edge, and critical curling stress occurs at the interior. 


e For both Westergaard and sensitivity analysis, the critical loading stress appears at the 
edge, and critical curling stress occurs at the interior. 
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Computed stresses with Westergaard method are significantly lower than KENPAVE 
software stresses for traffic loading whereas the difference for curling stresses is less. 


Sensitivity analysis is a significant tool in finding the optimum design of rigid pavement 
structure. 


It is recommended that the optimum thickness of the slab is 300 mm for selected values 
of input parameters as loading and curling stresses are minimum at this depth. 


KENPAVE software program provides more accurate results for stress analysis as 
Westergaard method has certain limitation while considering underneath support to the 
concrete slab. 
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